Does a Fermi liquid on a half-filled Landau level have Pomeranchuk instabilities? 
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We present a theory of spontaneous Fermi surface deformations for half-filled Landau levels (fill- 
ing factors of the form f = 2n + 1/2). We assume the half-filled level to be in a compressible, 
Fermi liquid state with a circular Fermi surface. The Landau level projection is incorporated via a 
modified effective electron-electron interaction and the resulting band structure is described within 
the Hartree-Fock approximation. We regulate the infrared divergences in the theory and probe the 
intrinsic tendency of the Fermi surface to deform through Pomeranchuk instabilities. We find that 
the corresponding susceptibility never diverges, though the system is asymptotically unstable in the 
n oo limit. 
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The discoveries of the integer quantum Hall effect 
(IQHE) 111 and the fractional quantum Hall effect 
(FQHE) y have stimulated many studies on the prop- 
erties of two-dimensional (2D) strongly correlated elec- 
tronic systems in a magnetic field. New ideas such as 
existence of fractionally charged quasiparticles 0] or the 
prediction of composite fermions '4*] have broadened our 
understanding of nature and have profoundly affected 
physics and other sciences. Theories built on these ideas 
have also proven very reliable to explain unusual prop- 
erties of strongly correlated electronic systems in the ex- 
treme quantum Hall regime at various filling factors. 

Somewhat more poorly-understood are the half-filled 
states, in particular the anisotropic states in high Lan- 
dau level-s (LL-s) 5]. The issue of anisotropy was ad- 
dressed theoretically quite early on [1, 0]- Moessner 
and Chalker 0] showed using Hartree-Fock theory that a 
striped charge density wave (CDW) prevails in the limit 
of very high LL-s if the interaction is of sufficiently short 
range relative to the magnetic length. In fact ample theo- 
retical evidence of their existence has accumulated ^ Q . 

While a stripe CDW phase is sufficient to explain 
anisotropy, other, competing structures may lead to an 
anisotropic response as well. Indeed it was pointed out 
in Ref. IT'I that, when the range of the effective interac- 
tion is comparable to the magnetic length, uniform states 
may win over the striped phase and, a few years later, it 
was shown [T^] that melting of the stripes could lead to 
a nematic phase. In this state, translational symmetry is 
completely restored but the system remains anisotropic. 
Interestingly such smectic and nematic states can be re- 
garded as the 'missing links' between the Wigner crystal 
and Fermi liquid states, based on a proposed general pic- 
ture of strong correlations (for a brief overview, see [ll|). 

In order to explain the emergence of anisotropy in 
Fermi liquid states at filling factor = 2n+ 1/2 (n > 2), 
there has been a surge of interest in the Pomeranchuk 



instability (PI) [l2|. Through such mechanism, a com- 
pressible Fermi liquid state (presumably the half-filled 
states in high LL-s) may " spontaneosly" enter into an 
anisotropic nematic state characterized by a deformed 
Fermi surface. In this scenario anisotropy emerges at one- 
particle level. In fact, the wave function for the nematic 
state proposed in Ref. 13] consists of single-particle 2D 
plane- wave states that form an elliptical Fermi sea 14 1. 
It is also plausible to suggest that anisotropy may emerge 
at a two-particle level, too. In this second scenario, 
one would start with a broken rotational symmetry wave 
function that contains a suitable symmetry-breaking pa- 
rameter in the two-particle correlation part of the wave 
function and a Slater determinant of 2D plane waves that 
form a standard circular Fermi sea 

In this Letter we test the first of the above two scenar- 
ios. Specifically, we address the question of formation of a 
nematic state (or its higher angular momentum generali- 
sations) from the opposite direction in the phase diagram 
of Ref. i.e. from the Fermi liquid side. Assuming a 
circular Fermi liquid ground state as a starting point, our 
modus operandi is to derive an effective electron-electron 
interaction that takes into account the projection onto 
the relevant half-filled LL. We then investigate whether 
this circular Fermi surface is unstable to small, point 
group symmetry-breaking deformations that are the 2D 
counterpart of the three-dimensional (3D) PI [l3|. 

As was discussed more generally in 3D, for several 
model effective interactions 17|, key requirements to 
have a PI are: (i) a sharp feature in the interaction po- 
tential at some characteristic distance; and (ii) for this 
characteristic length to be larger than the separation 
distance between particles in the system. Indeed the 
range of the effective interaction potential that we de- 
rive is longer than the magnetic length. Moreover it has 
non-monotonic features that develop into a sharp kink 
at a particular distance that increases as the order of 
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the LL-s increases. On the other hand, the average dis- 
tance between particles in a half-filled LL is fixed by the 
mangetic length. Thus both of the above conditions are 
met asymptotically for sufficiently high LL-s. Testing the 
ensuing expectation of an intrinsic tendency to a PI for 
high half-filled LL is the main purpose of this work. 

We consider filling factors of the form: v = 2n + v* 
where n = 1,2,... is the index of the uppermost half- 
filled LL {v* = 1/2) assumed to be fully spin-polarized. 
According to the Halperin-Lee-Read theory [l8|, the N* 
electrons of the half-filled LL form a 2D circular Fermi 
liquid state (with density p*) that effectively sees no 
magnetic field. We have i^* — 2'kIqP* where the mag- 
netic length /q = \/hc/eB represents the characteris- 
tic length scale in the problem. Since v* = 1/2 and 
p* — TT fc|,/(2 tt)^, the radius of the circular 2D Fermi 
wave vector is kp = I/Iq. Assuming that kinetic energy 
is quenched the Hamiltonian, up to a constant, is 

Hn^lJ d\ J dVV„(r)i;(|f-f'|)p„(r") , (1) 

where v{\f — r'|) = e^/(4 7re|r — r'\) is the Coulomb 
potential. For simplicity we will take the dielectric func- 
tion e to be a constant. Here p„(r) = \E'}i(r)5'„(r) is the 
density operator. The quantum field operator, 'i'l^{f), 
creates a spinless electron in the n^'' LL at position r. 

It is well known that projection onto the n}^ LL con- 
tained implicitely in the definition of the density oper- 
ator Pnif) introduces highly non-trivial physics in the 
Hamiltonian in Eq. ([1]). The key idea of our approach 
is to assume a uniform circular Fermi liquid state and 
investigate possible instabilities induced by the heavily 
renormalised interactions. 

One can transform the density operator in 2D Fourier 
space as Pn{q) = J <Pr e^^^ pn{r) and rewrite Eq.lH]) as 
Hn = ^J2q'"il) Pni-q) Pn{q} whcrc v{q) = 2Tre'^ / {Aneq) 
is the 2D Fourier transform of the Coulomb interaction 
potential. Based on the Hamiltonian theory of Ref. [l^ 
we can project the density onto the half-filled LL, which 
amounts to writing Pn{q) — Fn{q)'p{q} where ]}{(]) is the 
projected density operator. Thus the Hamiltonian can 
be rewritten in 2D Fourier space as 



(2) 



where Vniq) = Fn{q)'^ v{q) represents an effective inter- 
action potential which in real space would be given by 
Vn{r) = / (Pq/{2TTf e-^^^Vn{q). The form factors are 
Fn{q) — Ln (9^^0/2) exp (— 5^/q/4) where the L„(a;)-s are 
Laguerre polynomials. 

There are two consequences of the projection onto the 
half-filled Landau level implicit in the Hamiltonian ([2]): 
firstly, the effective interaction potential Vn(r) is heav- 
ily renormalized when compared to the Coulomb inter- 
action. Secondly, the projected density operators them- 



selves have a non-trivial algebra. Here we wish to in- 
vestigate whether the first of these two features may be 
sufficient to produce a PI instability in a putative Fermi 
liquid state. In line with this view, we assume that new 
fermion creation and annihilation operators "0^1 "0 can 
be introduced in such a way that the projected densi- 
ties can be written as p{q) — V'(^)^0(^ + (a form 
that implies standard operator algebra). Here ip^k) cre- 
ates a fermion in a plane wave state. One can then ap- 
ply the mean-field theory of a PI in a 3D continuum 
in Ref. [l3| (recently generahsed to 2D [13, [13). It 



starts with a Hartree-Fock ansatz for the ground state, 

\h) = rife ^h^^'^ io) 

the vacuum. Note that this ansatz is a homogeneous, 
itinerant state with an arbitrary dispersion relation e^. 
For the Hamiltonian in this Slater determinant of 
plane waves affords a rudimentary description of the re- 
emergence of itinerancy when the kinetic ener gy h as been 
completely quenched (note that, unlike Refs. 17, 2^, 21|, 
there is no 'bare' contribution to e^r here - it all comes 
from interactions) . The functional form of is our varia- 
tional parameter. It determines which plane wave states 
are occupied via < and is found by minimization 
of {Hn)- That yields a self-consistency equation giving 
the shape of the Fermi surface through e^^ = 0. PI-s are 
point group symmetry-breaking instabilities of the shape 
of the Fermi surface. 

To find an instability equation, it is natural to split 
in two parts: one that preserves the continuous rotational 
symmetry of the plane and another one that may break 
it. We thus write ill 



^eo{\k\)-ki{\k\)cos{ie^), 



(3) 



where eod^l) is the symmetric component of the disper- 
sion relation and ^ = 1, 2, 3, . . . determines the symmetry 
of the instability. The condition of instability towards a 
small deformation of the Fermi surface is (ignoring the 
possibility of a first-order phase transition) [2O, [2l| 



(4) 



where V/ = At: I dr rVnir) Ji{k"p rY (5) 

measures the strength of the electron-electron interaction 
in the channel with angular momentum Z = 1, 2, 3, . . . For 
I ~ I, the PI corresponds to a rigid displacement of the 
Fermi suface in reciprocal space, without change of ei- 
ther shape or volume. This can never lead to a lowering 
of the energy in a Galilean-invariant system and there- 
fore this instability cannot take place [l3, [Hi . [This 
also follows from the more explicit form of the instability 
equation derived below, Eq.((7]).] On the other hand for 
Z = 2, 3, 4 ... we could have PTs corresponding to defor- 
mations of the Fermi surface possessing d-wave, /-wave. 
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FIG. 1: Interaction potentials calculated from Eq.Q along- 
side the large-n [Eq. ([9}] and large-^ (Coulomb) asymptotes. 



g-wave... symmetry, respectively. The Fermi velocity in 
Eq. (gl) is given by 
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FIG. 2: Dependence of I(n, I) in Eq. (fTTjl on the cutoff ("c 
for w — 1 (red solid curves), 2 (green long dash) and 4 (blue 
short dash) for (a) n = 1; (b) n = 2; and (c) n = 10 (note the 
different scales on the third plot). We employed the asymp- 
totic expression in Eq. © to plot (c). The values of angular 
momentum I are as indicated. 
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Vi, 
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where we take into account the infinite bare mass m oo 
of our Hamiltonian [see Eq. Substituting this into 

the instabiUty equation (U) yields 



Vi-Vi>0. 



The effective interaction can be written as 



(7) 



Vn{r) 
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(8) 

where Ji{x) are Bessel functions and n = 1, 2, 3, . . . We 
can write Vn{r) — {e^ / Airelo) Vn{C) where ( — r/lo is the 
natural dimensionless distance. We calculated fn(C) ex- 
actly for several n (we do not display such expressions 
for the sake of brevity). The results are shown in Fig. [U 
where we also plot an asymptotic expression for the ef- 
fective interaction obtained in the limit of high LL-s: \^ 
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n > 1. (9) 



K{x) is the complete elliptic integral fo the first kind. All 
interaction potentials tend symptotically to the Coulomb 
interaction, ~ 1/C for large C,. This leads to in- 

frared divergences in the theory whose regulation is dis- 
cussed below. On the other hand, the presence of a length 
scale in the problem (the magnetic length) is apparent at 
shorter distances where one notices that as n increases an 
increasingly sharp kink develops at the specific distance 
Tn = Cnlo where 



C„ = 2V2n 



(10) 



Note in particular that the ultraviolet divergence of the 
Coulomb interaction at short distance has been sup- 
pressed. Interestingly, a sharp feature in the interaction 
potential at a finite distance r„ suggests the possibility 
of a PI provided the dimensionless parameter rnkp is 
large enough [17|. Since Vnkp — (n, Eq. pO)l imphes an 



increased tendency towards a PI in high LL-s. 

Note that since Vn (C) 1/C for C > 1 the RHS of 
Eq. (O diverges. This leads to diverging Landau param- 
eters, Fi cx Vi [2l| and Fermi velocity Vp oc Vi. These 
infrared divergences can be regulated by introducing a 
cutoff Cc in the electron-electron interaction potential. 
The actual values of vp and Fi depend on the actual 
value and form of the cutoff which in turn depends on 
extrinsic features of the system such as sample size or 
device configuration. However, as long as the cutoff dis- 
tance is long enough, the instability condition does not 
depend on these parameters. It is such intrinsic ten- 
dency to a PI that we wish to probe further. To do so let 
us rewrite the instability condition in Eq.([7]) in a more 
explicit form: 

i{n,i) ^ / dCC«n(C) [MO' - MO'] > . (11) 

Jo 

Clearly, if /(n, Z) < 0, a PI does not occur. However, 
as I{n, I) becomes less negative, the susceptibility to a 
PI with angular momentum quantum number I increases 
until it diverges at /(n, I) = 0. The following replacement 



Vn (C) - / ( ) Vn (C) , 



(12) 



with / (x) = jFqrj- provides a mathematically convenient 
way of introducing a cutoff in the form of a smoothed-out 
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FIG. 3: Dependence of the integral I(n,l) in Eq. eval- 
uated numerically, on n. The angular momenta of the insta- 
bilities are / = 2, 3, 4 and 5 (solid, long-dashed, short-dashed 
and dotted lines, respectively). Circles represent values ob- 
tained using the effective potentials, i'n(C) for n = 1,2 and 
3. Diamonds were obtained using the asymptotic formula in 
Eq.©. The inset shows the behaviour at very large n for the 
I = 2 case. The straight line in the inset is a fit of the form 
I{n,l) — Ai/n^' to the numerical data for n > 800. The 
fitting parameters are A2 « 0.15 and B2 ~ 0.54. 

step function of finite width w. Fig. [5] shows the depen- 
dence of I{n, I) on the cutoff (c- For sufficiently large 
the integral is independent of Q as well as of the width 
w of the cutoff. Taking the limit Cc ^ 00 represents 
a convenient way to estimate the intrinsic tendency of 
the system to a PI which is independent of the form and 
value of the cutoff. Physically, such a cutoff may corre- 
spond to, for example, the finite thickness of the device. 
The crucial point is that the value of the integral I{n, I) 
in Eq. pip is an intrinsic feature of the system inde- 
pendent of the form, size and mechanism of the cutoff 
as long as is much larger than the average separation 
between particles ~ kp^ = Iq, and ^ Cn where C„ is 
the distance where the kink of the effective potential in 
Fig. [T] appears. The dependence of the converged values 
of I{n, I) on n and / is shown in Fig. [3] We note that 
the PI condition is never met for v = 5/2, 9/2 and 13/2. 
The susceptibility to a PI has a non-trivial dependence 
on n with a maximum for I = 2 at n — 2, while it mono- 
tonically increases for Z > 2 up to n = 3. Also note that 
the system seems most susceptible to a PI in the I — 2 
channel. From these results we conclude that the system 
never has an intrinsic PI (i.e. /(n, 2) is always < 0). On 
the other hand, for large n we have I{n,2) (x 1/n'^-^'*, 
which implies that the half-filled level is asymptotically 
unstable to a PI in the n ^ 00 limit [2^ . 

In summary, we have investigated whether a compress- 
ible circular Fermi liquid state in half-filled high LL-s un- 
dergoes a phase transition to a non-circular anisotropic 
nematic phase through a PI. We use the Hamiltonian 



theory approach [l^ to derive a properly LL projected 
effective electron-electron interaction from which itiner- 
ancy emerges at the mean-field level. We look for defor- 
mations of the Fermi surface by testing for an intrinsic 
divergence of the corresponding susceptibility. We find 
that the susceptibility towards a PI is increasingly large 
as we move to high LL-s and diverges in the n — > cxd 
limit. The increased tendency towards a PI is a direct 
consequence of the length scale r„ ^ ^/2nlo present in 
the effective interaction. 

As is well known, there is also a tendency towards 
stripe formation 0] with whom the PI competes. In 
real systems, there is a second length scale, Tc, responsi- 
ble for cutting off infrared diveregences associated with 
the long-range nature of the Coulomb interaction. Ev- 
idently r„ and r^ would become comparable for large 
values of n. The intrinsic susceptibility to a PI is already 
quite high in this regime. Thus, even minor extrinsic 
effects (sample thickness, device configuration, etc.) as- 
sociated with this second length scale may trigger such 
an instability. A calculation of the critical value of n 
at which this transition might take place is beyond the 
scope of the present analysis. It would require taking 
device configuration into account as well as comparing 
the energy of any Fermi liquid state to those obtained 
for stripe configurations. 
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